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The influence of the intramolecular anharmonicity and the 
strong vibron-phonon coupling on the two-vibron dynamics 
in an a-helix protein is studied within a modified Davydov 
model. The intramolecular anharmonicity of each amide-I vi- 
bration is considered and the vibron dynamics is described 
according to the small polaron approach. A unitary trans- 
formation is performed to remove the intramolecular anhar- 
monicity and a modified Lang-Firsov transformation is ap- 
plied to renormalize the vibron-phonon interaction. Then, a 
mean field procedure is realized to obtain the dressed anhar- 
monic vibron Hamiltonian. It is shown that the anharmonic- 
ity modifies the vibron-phonon interaction which results in 
an enhancement of the dressing effect. In addition, both the 
anharmonicity and the dressing favor the occurrence of two 
different bound states which the properties strongly depend 
on the interplay between the anharmonicity and the dressing. 
Such a dependence was summarized in a phase diagram which 
characterizes the number and the nature of the bound states 
as a function of the relevant parameters of the problem. For 
a significant anharmonicity, the low frequency bound states 
describe two vibrons trapped onto the same amide-I vibration 
whereas the high frequency bound states refer to the trapping 
of the two vibrons onto nearest neighbor amide-I vibrations. 
PACS: 03.65.Ge, 63.20.Ry, 63.22.+m, 87.15.-v 



I. INTRODUCTION 

In low-dimensional molecular lattices, the nonlinear 
nature of vibrational excitons (vibrons) plays a key role 
for energy transfer as well as energy storage in both phys- 
ical, chemical and biological systems. In a general way, 
two main sources yield a nonlinear dynamics, namely the 
intrinsic intramolecular anharmonicity of each molecule 
and the extrinsic coupling between the vibrons and sur- 
rounding low frequency excitations (such as phonons for 
instance). 

In classical lattices, the anharmonicity gives rise to 
the occurences of intrinsic localized modes, or discrete 
breathers, which have been the subject of intense re- 
searches during the last decade (for a recent review see for 
instance Refs. [1-3]). These highly localized vibrations 
do not require integrability for their existence and cor- 
respond to quite general and robust solutions [4], Since 
discrete breathers favor a local accumulation of the en- 
ergy which might be pinned in the lattice or may travel 
through it, they are expected to be of fundamental im- 
portance for both energy storage and transport. Unfor- 



tunately, no clear evidence has yet been found for the 
existence of breathers in real molecular lattices. By con- 
trast, in the quantum regime, two-vibron bound states 
(TVBS) have been observed in several low-dimensional 
systems [5-13]. In that case, the intramolecular anhar- 
monicity breaks the vibron independence and favors the 
formation of bound states [14-19]. When two vibrons 
are excited, a bound state corresponds to the trapping of 
two quanta over only a few neighboring molecules with 
a resulting energy which is lesser than the energy of two 
quanta lying far apart. The lateral interaction yields a 
motion of such a state from one molecule to another, thus 
leading to the occurrence of a delocalized wave packet 
with a well-defined momentum. As a result, TVBS are 
the first quantum states which experience the nonlincar- 
ity and can thus be viewed as the quantum counterpart 
of breathers or soliton excitations [16]. 

The second source of nonlinearity, which originates in 
the coupling between vibrons and low frequency excita- 
tions, was first pointed out by Davydov and co-workers 
[20] to explain the mechanisms for bioenergy transport. 
The main idea is that the energy released by the hy- 
drolysis of ATP can be stored in the high frequency 
C=0 vibration (amide-I) of a peptide group of a pro- 
tein. The dipole-dipole coupling between the different 
peptide groups leads to the derealization of the internal 
vibrations and to the formation of vibrons. However, the 
coupling between the vibrons and the phonons of the pro- 
tein yields a nonlinear dynamics which counterbalances 
the dispersion created by the dipole-dipole interaction. 
It leads to the creation of the so called Davydov's soli- 
ton which provides an approximation to the self-trapping 
phenomena described by a Frohlich type Hamiltonian 
[21]. Soliton mechanisms for bioenergy transfer in pro- 
teins have received increasing attention during the last 
twenty five years and a broad review can be found in 
Refs. [22,23]. 

However, as pointed out by Brown and co-workers 
[24,25] and by Ivic and co-workers [26-28], the solution of 
the Davydov's problem is rather a small vibron-polaron 
than a vibron-soliton. Indeed, the self-trapping process 
exhibits two asymptotic solutions depending on the val- 
ues taken by the three relevant parameters of the prob- 
lem, i.e. the vibron bandwidth, the phonon cutoff fre- 
quency and the small polaron binding energy propor- 
tional to the strength of the vibron-phonon coupling. 
When the vibron bandwidth is greater than the phonon 
cutoff frequency, the adiabatic limit is reached. The 
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phonons behave in a classical way and create a quasi- 
static potential well responsible for the trapping of the 
vibron. The vibron, dressed by this lattice distortion, 
forms a polaron which is described according to the soli- 
ton theory of Davydov. By contrast, when the vibron 
bandwidth is lesser than the phonon cutoff frequency, the 
situation corresponds to the non-adiabatic limit in which 
the quantum behavior of the phonons plays a crucial role. 
As mentioned by the authors, this situation corresponds 
to the vibron dynamics in proteins. Therefore, a vibron 
is dressed by a virtual cloud of phonons which yields a 
lattice distortion essentially located on a single site and 
which follows instantaneously the vibron. The dressing 
effect modifies the vibron frequency, reduces the vibron 
bandwidth and allows for an attractive interaction be- 
tween vibrons mediated by virtual phonons. Such an in- 
teraction is responsible for the formation of bound states 
and it has been suggested that proteins can support soli- 
tons formed by bound states involving a large number of 
vibrational quanta [27,28] . 

As a consequence, the previous results clearly show 
that both the intramolecular anharmonicity and the 
strong vibron-phonon coupling produce a similar effect 
onto the vibron dynamics and favor the formation of 
bound states. The present paper is thus devoted to the 
fundamental question of the interplay between both non- 
linear sources. To proceed, the Davydov model, mod- 
ified by introducing the intramolecular anharmonicity 
of each amide-I vibration, is described within the small 
polaron approach. The two-vibron energy spectrum is 
studied with a special emphasis on the influence of the 
intramolecular anharmonicity onto the dressing effect. 
Note that we do not investigate the formation of solitons 
but focus how attention onto the creation of TVBS, only. 
Such a procedure is twofold. First, as mentioned previ- 
ously, TVBS are the first quantum states sensitive to the 
nonlinearity. Their characterization allows us to under- 
stand the interplay between both nonlinear sources and 
appears as a first step to study the formation of multi- 
vibron solitons. Then, recent theoretical calculations 
have shown that single-vibron solitons do not last long 
enough to be useful at biological temperatures [29,30]. 
By contrast, two-vibron solitons are more stable and ap- 
pear as good candidates for bioenergy transport [31,32]. 
Note that a perfect knowledge of the two-vibron dynam- 
ics is also required to interpret some experiments such as 
time resolved pump-probe spectroscopy [33,34]. 

The paper is organized as follows. In Sec. II, the Davy- 
dov Hamiltonian for a one-dimensional molecular lattice 
is described. In Sec. Ill, we first realize a unitary trans- 
formation to remove the intramolecular anharmonicity. 
Then, a modified Lang-Firsov [35] transformation is ap- 
plied to renormalize the vibron-phonon interaction and 
to reach the small polaron point of view. Finally, a mean 
field procedure is performed to obtain the dressed anhar- 
monic vibron Hamiltonian. In Sec. IV, we summerizc 
the number states method used to solve the two-vibron 
Schrodinger equation. In Sec. V, a detailed analysis of 



the two-vibron energy spectrum is performed depending 
on the values taken by the relevant parameters of the 
problem. Finally, the results are discussed and inter- 
preted in Sec. VI. 



II. THE DAVYDOV HAMILTONIAN 

According to the original Davydov model, the vibron- 
phonon dynamics in an a-helix protein is described by 
a one-dimensional system formed by N sites periodically 
distributed along the lattice. Each site is occupied by 
a peptide group which contains the amide-I vibration. 
The nth amide-I vibration is assumed to behave as an 
internal high frequency oscillator described by the stan- 
dard creation and annihilation vibron operators 6+ and 
b n . The vibron Hamiltonian is thus written as (using the 
convention h=l) 

H V = J2 ^ob+b n + 73 (6+ + b n f + 74(6+ + b n f 

n 

-J(b+ + b n )(b+ +1 + b n+1 ) (l) 

where wo stands for the internal frequency of the nth 
amide-I mode and where J denotes the lateral hopping 
constant between nearest neighbor amide-I vibrations. In 
Eq.(l), 73 and 74 represent the cubic and quartic anhar- 
monic parameters of each amide-I mode. 

The amide-I vibrations interact with the phonons of 
the lattice which characterize the collective dynamics of 
the external motions of the peptide groups. Within the 
harmonic approximation, each peptide group, with mass 
M, interacts with its nearest neighbor peptide groups 
via the lateral force constant W. Therefore, the phonons 
correspond to a set of N low frequency acoustic modes, 
labeled {q} 7 which the Hamiltonian is defined as 



Hp = ^2 £l q a q a q 



(2) 



where a q and a q stand for the phonon operators of the 
qth. mode with frequency Q q = Cl c \ sin(q/2) |, Q c = 
y^AW/M denoting the phonon cutoff frequency. 

Finally, the vibron-phonon interaction Hamiltonian, 
which characterizes a random modulation of the inter- 
nal frequency of each amidc-I mode, is expressed as 



Aif„ 



1 



^ 2 

qn 
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where the coupling constant A„ i9 is written as 
. Aq sin(q) 



'VN y/\ sin(q/2) 



-iqn 



(3) 



(4) 



In Eq.(4), Ao is defined in terms of the coupling pa- 
rameter x introduced in the original Davydov model as 
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A = x(?i 2 -MW) -1 / 4 (h as be reintroduced to avoid con- 
fusion) . 

The vibron-phonon dynamics is thus described by the 
full Hamiltonian H = H v + H p + AH vp which slightly 
differs from the original Davydov model. The Davydov 
Hamiltonian is recovered from Eqs.(l)-(3) by restricting 
the full Hamiltonian to vibron-conserving terms, only, 
and by neglecting the intramolecular anharmonicity of 
each amide-I vibration. Although the full Hamiltonian 
H yields a rather simple model for the protein dynamics, 
it cannot be solved exactly due to the anharmonic con- 
tributions. The following section is thus devoted to its 
simplification to obtain an effective Hamiltonian describ- 
ing the vibron dynamics. 

III. DRESSED ANHARMONIC VIBRONS 

A. Renormalization of the intramolecular 
anharmonicity 

To remove the cubic and quartic intramolecular anhar- 
monicity, the standard procedure described by Kimball 
et al. [14] is used. First, by disregarding the lateral cou- 
pling between nearest neighbor amide-I vibrations as well 
as the vibron-phonon interaction, the procedure consists 
in performing a unitary transformation T to diagonalize 
each anharmonic amide-I mode. Then, an approximate 
Hamiltonian is obtained by applying the transformation 
on the full Hamiltonian H and be keeping the vibron- 
conserving terms, only. 

As a result, by using the formula shown in appendix A, 
the transformed vibron-phonon Hamiltonian H = THT + 
is expressed as 

H = ]T(co -2A- B)b+b n - Ab+ 2 b 2 n - Bb+ +1 b+b n+1 b n 

n 

- Ji[6+6 n+1 + h.o.] - J 2 [b+ 2 b 2 n+l + h.o] 

- M b n( b n b n + b+ +1 b n+1 )b n+1 + h.o] 

+ Y,( A ^ a t + Kq^K 1 + 2vKb n + V b+ 2 b 2 n ] 

qn 

+ J2 n i a t a i ( 5 ) 

where A denotes the positive anharmonic parameter writ- 
ten as 

A = 30— 2 - 674 (6) 

UJ 

In Eq.(5), h.o. stands for the corresponding hcrmitian 
operator and the different parameters are defined as 

B = 144J(^) 2 ; Ji = J(l + 44(^) 2 - 12^) 

J2 = 4J(— ) 2 ; J 3 - J(22(^) 2 - 12^) 

v = 120(^) 2 - 12^* (7) 
ujq ujo 



Note that Eq.(5) was obtained by disregarding the con- 
stant term as well as contributions acting onto states in- 
volving more than two vibrons. 

The unitary transformation allows us to diagonalize 
each anharmonic amidc-I vibration up to the second or- 
der with respect to the anharmonic parameters. As a re- 
sult, the vibron operators describe new vibrational states, 
called anharmonic vibrons, which are expressed as linear 
superimpositions of the harmonic states of each amidc- 
I mode. Therefore, Eq.(5) clearly shows that the an- 
harmonicity strongly affects the dynamics of these an- 
harmonic vibrons. It first modifies the harmonic part 
of the Hamiltonian by inducing a redshift of each inter- 
nal frequency (u>o — — B) and by changing the 
strength of the lateral interaction (J — > J\). Then, the 
anharmonicity is responsible for the occurrence of cou- 
pling terms which break the vibron independence and 
which directly affect the two-vibron dynamics. The terms 
—Ab+b+b n bn and — Bb^b^ +1 b n b n +i yield an attractive 
interaction between two vibrons and favor their trapping 
around the same amide-I site and around two nearest 
neighbor sites, respectively. The contributions propor- 
tional to J2 characterize hops in the course of which two 
vibrons, initially located on the same amide-I vibration, 
realize simultaneously a transition to a nearest neighbor 
amide-I site. In the same way, the contributions propor- 
tional to J3 affect single- vibron hops from states formed 
by two vibrons located onto the same amide-I vibration. 
Finally, the anharmonicity leads to a correction of the 
vibron-phonon coupling Hamiltonian and Eq.(5) shows 
that single- and two-vibron states do not experience the 
same interaction with the phonon bath. When com- 
pared with the harmonic situation, the coupling between 
phonons and singlc-vibron states is enhanced by the an- 
harmonicity, i.e. it is multiplied by the factor 1 + 2ry. 
Moreover, the terms r/b^b^b n b n , which act onto two- 
vibron states only, show that two vibrons located onto 
the same amidc-I mode interact in a different way with 
the phonons bath when compared with two vibrons which 
are far apart. Note that the anharmonic corrections oc- 
curring in Eq.(7) are typically of the order of A/u>o and 
represent rather small contributions which have been ne- 
glected in our previous works [17-19]. Nevertheless, al- 
though some changes can effectively be disregarded, such 
as the modification of the hopping constants, other con- 
tributions play a significant role, especially the correction 
of the vibron-phonon interaction, as it will be shown in 
the following sections. 

B. Renormalization of the vibron-phonon coupling 

The next step in our procedure is to partially remove 
the vibron-phonon coupling Hamiltonian by performing 
a modified Lang-Firsov [35] transformation. According 
to the Ivic and co-workers's remarks [27,28], the vibron- 
phonon dynamics is dominated by the so-called dress- 
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ing effect since the non-adiabatic limit is reached. As a 
result, we consider a "full dressing" and introduce the 
following unitary transformation 

U = exp(Y, Qn[(l + 2 V )b+b n + V b+b+b n b n }) (8) 

n 

where Q n is defined as 

g» = D^tf - tI^] (9) 

The transformation U (Eq.(8)) slightly differs from the 
true Lang-Firsov transformation due to its dependence 
with respect to the vibron population operators. This 
dependence originates in the modification of the vibron- 
phonon interaction mediated by the intramolecular an- 
harmonicity (see Eq.(5)) and the original Lang-Firsov 
transformation is recovered by setting 77 = 0. There- 
fore, by using Eq.(8), the transformed Hamiltonian H = 
UHU + is expressed as 

H = J2^ob+b n - Ab+ 2 b 2 n - Bb+ +1 b+b n+1 b n 

n 

- Ji[e+(7v„ - i)e n+ i(jv n+ i)6+& n+ i + h.o] 

- J 2 [Q + n 2 (N n \)Q 2 n+l (N n+1 + \)b+%l +l + h.o.] 

- M®n(N n - l)Qn+l(N n+1 )b+[N n + N n+1 ]b n+1 + h.o] 

+ J2n q a+a q (10) 

where N n = b^b n denotes the vibron population operator 
and where the different parameters occurring in Eq.(10) 
are expressed in terms of the small polaron binding en- 
ergy E B = 2Al/Q c as 

uj = uj - 2A - B - (1 + Ar])E B 
A = A+(1 + 8t])E b ;B = B + (1 + 4 V )E b (11) 

In Eq.(10), 6„(iV n ) stands for the dressing operator de- 
fined as 

6„(JV„) = exp{-Q n [l + 2r ] + 2r 1 N n }) (12) 

In this dressed anharmonic vibron point of view 
(Eq.(10)), the vibron-phonon coupling remains through 
the modulation of the lateral terms by the dressing oper- 
ators ® n (N n ). Although these operators depend on the 
phonon coordinates in a highly nonlinear way, the vibron- 
phonon interaction has been strongly reduced within this 
transformation. As a result, we can take advantage of 
such a reduction to express the full Hamiltonian H as 
the sum of three separated contributions as [26,27] 

H = H e ff + H p + AH (13) 

where H e ff =< H > — H p denotes the effective Hamil- 
tonian of the dressed anharmonic vibrons and where 



AH = H— < H > stands for the remaining part of the 
vibron-phonon interaction. The symbol < ... > repre- 
sents a thermal average over the phonon degrees of free- 
dom which are assumed to be in thermal equilibrium at 
temperature T. As a result, the effective dressed anhar- 
monic vibron Hamiltonian is written as 

H eff = J2^b n - Ab+ 2 b 2 n - Bb+ +1 b+b n+1 b n 

n 

- Ji[$(JV„ + N n+1 )b+b n+1 + h.o.} 

-JiMN n + N n+1 ) 4 b+ 2 b 2 n+1 + h.o.} 

-J 3 [$(N n + N n+1 )b+[N n + N n+1 }b n+1 + h.o.} (14) 

where $(A) = exp(-S(T)[l + 2r] + 2r/X}) and where 
S(T) is the coupling constant introduced by Ivic and co- 
workers as (k B denotes the Boltzmann constant) 

q 

The Hamiltonian H e ff (Eq.(14)) describes the dynam- 
ics of the anharmonic vibrons dressed by a virtual cloud 
of phonons, i.e. anharmonic small polarons. It takes 
into account on the anharmonicity up to the second or- 
der and allows for a renormalization of the main part 
of the vibron-phonon coupling within the non-adiabatic 
limit. The interaction Hamiltonian AH, which charac- 
terizes the coupling between these dressed anharmonic 
vibrons and the remaining phonons, is assumed to be 
small in order to be treated using perturbation theory. 
Such a contribution, responsible for phase relaxation, will 
be studied in a forthcoming paper. 

Eq.(14) clearly shows the interplay between the in- 
tramolecular anharmonicity and the strong vibron- 
phonon coupling. The Hamiltonian H e ff exhibits, ba- 
sically, the same contributions as the vibron part of the 
anharmonic Hamiltonian H. The main difference is that 
the parameters occurring in H e f f are renormalized due to 
the dressing effect which modifies the dynamics accord- 
ing to two main ways. First, it yields additional contribu- 
tions to the anharmonic parameters A and B as well as to 
the internal frequency u>o (Eq.(ll)). Then, the dressing 
effect modifies the different lateral contributions via the 
dressing function &(X) which reduces the hopping con- 
stants. However, in a marker contrast with the harmonic 
situation, the intramolecular anharmonicity enhances the 
role played by both the small polaron binding energy and 
the coupling constant S(T). In addition, the dressing ef- 
fect depends on the vibron population. As shown in the 
next section, such a dependence discriminates transitions 
involving two vibrons located onto the same amide-I vi- 
bration from transitions involving two vibrons which are 
far apart. 
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IV. TWO- VIBRON STATES 

Since the effective Hamiltonian H e ff (Eq.(14)) con- 
serves the number of dressed anharmonic vibrons, its 
eigenstates can be characterized by using the number 
states method [16-19]. Within this method, the two- 
vibron wave function is expanded as 

I *> = H *(ni,«2) I n u n 2 ) (16) 

ni ,n2 >7li 

where {| m, n 2 )} denotes a local basis set normalized and 
symmetrized according to the restricting n 2 > n\ due to 
the indistinguishable nature of the vibrons [17-19]. A 
particular vector | 711,712) characterizes two vibrons lo- 
cated onto the sites 71 1 and n 2 , respectively. The expres- 
sion of the corresponding time independent Schrodingcr 
equation, H e ff I = w I depends on the nature 
of the basis vectors involved in. Indeed, from Eq.(14), 
there are three different situations either describing two 
vibrons located onto sites which are far apart, two vi- 
brons located onto nearest neighbor sites and two vibrons 
located onto the same site. 

In the first situation, the Schrodinger equation is ex- 
pressed as 

- Ji$(l)[*(ni,n 2 + 1) + *(ni,n 2 - 1)] 

- Ji$(l)[*(ni + 1, n 2 ) + *(ni - l,n 2 )] 

+ 2a> *(ni,n 2 ) =w*(ni,n 2 ) (17) 

Eq.(17) shows that the two vibrons move in an inde- 
pendent way according to an effective hopping constant 
Ji$(l) = Jtexp[-(1 + 4?7)S(T)] which slightly differs 
from the hopping constant involved in the harmonic ap- 
proximation. Indeed, as discussed in the previous section, 
the intramolecular anharmonicity modifies the harmonic 
hopping constant (J — > J\) and enhances the dressing 
effect (S(T) -> (1 + 4r])S(T)). Therefore, although the 
anharmonic parameter rj is rather small, it acts under 
the exponential and then reduces the effective hopping 
constant with respect to the harmonic situation. 

When two vibrons are located onto nearest neighbor 
amide-I vibrations, the Schrodinger equation is expressed 
as 

- Ji$(l)[*(m,ni + 2) + *(m - 1,77! + 1)] 

- V2(Ji + J 3 )*(2)[*(ni,m) + *(m + l,m + 1)] 

+ (2^o - B)*(ni,m + 1) = wtf (ni,m + 1) (18) 

and, when two vibrons are located onto the same amidc-I 
mode, the Schrodingcr equation is defined as 

- 2J 2 $(2) 4 [4'(ni + 1, 77! + 1) + 4/(771 - 1, 77i - 1)] 

- V2(Ji + J 3 )*(2)[*(ni,m + 1) + tf (m - l,m)] 

+ (2a>o-2i)*(ni,ni)=w*(m,ni) (19) 

As shown in Eqs.(18)-(19), both the intramolecular an- 
harmonicity and the dressing effect strongly modify 



the dynamics involving vibrons located onto neighbor- 
ing sites. First, they contribute to a redshift of the 
energies of the corresponding states. Then, they af- 
fect the hopping processes by favoring simultaneous mo- 
tions of two vibrons and by changing the dressing ef- 
fect which is characterized by the effective correction 
$(2) = exp[-(l+6r])S(T)]. This latter feature originates 
in the dependence of the dressing operators (Eq.(12)) on 
the vibron population. Therefore, since $(2) is lesser 
than $(1), it is clear that the trapping process experi- 
enced by two vibrons located onto the same site is more 
efficient than the dressing effect which modify the singlc- 
vibron dynamics. 

Eqs.(17)-(19) clearly indicate how the physics of the 
two-vibron states is related to the dynamics of a sin- 
gle fictitious particle moving quantum mechanically on 
the two-dimensional lattice displayed in Fig. la [17-19]. 
Within this equivalence, the two-vibron wave function 
4 , (t7i, 77 2 ) can be viewed as the wave function of the fic- 
titious particle. Its dynamics is described by a tight- 
binding Hamiltonian characterized by the self-energy 2u>o 
located on each site and a hopping matrix ^$(1) which 
couples nearest neighbor sites. However, the 2D lattice 
exhibits two rows of defects (see Fig. la) which yield a 
redshift of the self-energy of the corresponding sites as 
well as for a modification of the hopping matrix elements 
connecting the defect sites. Therefore, such defects allow 
us to discriminate between two different eigenstates. The 
eigenstates of the core of the lattice correspond to plane 
waves slightly perturbed by the defects. By contrast, the 
presence of the defects leads to the occurrence of states 
which are localized in the vicinity of the two defect rows. 
In terms of the two-vibron dynamics, the previous equiva- 
lence yields two different kinds of eigenstates. Indeed, the 
derealization of the fictitious particle is associated to a 
free motion of the two vibrons, i.e. two-vibron free states 
(TVFS), whereas its localization is connected to the oc- 
currence of TVBS. Note that in marked contrast with the 
standard Hubbard Hamiltonian for bosons [14,16,18], the 
presence of two defect rows in our equivalent lattice shows 
that the system can support two kinds of bound states 
corresponding to the trapping of the two vibrons either 
onto the same amide-I vibration or onto nearest neighbor 
amide-I sites, respectively. 

The Schrodinger equation Eqs.(17)-(19) can be ex- 
pressed in an improved way by taking advantage of the 
lattice periodicity. Indeed, the two-vibron wave function 
is invariant with respect to a translation along the lattice 
and can be expanded as a Bloch wave as 

*(ni, 7i 2 = 77! + 777) = -L V e ifc < ni+m / 2 >* fc (m) (20) 
\f N 

v ni 

In Eq.(20), the total momentum k, which belongs to the 
first Brillouin zone of the molecular lattice, is associated 
to the motion of the center of mass of the two vibrons 
whereas the resulting wave function \&fc(m) refers to the 
degree of freedom 777 which characterizes the distance 
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between the two vibrons. Since the momentum k is a 
good quantum number, the Hamiltonian H e ff appears 
as block diagonal and the Schrodinger equation can be 
solved for each k value. Therefore, when two vibrons are 
located onto sites which are far apart (m > 1), Eq.(17) 
becomes 

(2cj - w)**(m) - r fc [tf fc (m + 1) + * fe (m - 1)] = 

(21) 

where Tk = 2 Ji$(l)cos(fc/2). In the same way, when two 
vibrons are located onto nearest neighbor sites (m = 1), 
Eq.(18) is rewritten as 

(2lu -B- w)* fc (l) - r fe * fe (2) - V2 7fc * fc (0) = (22) 

where j k = 2(J 1 + J 3 )$(2)cos(fc/2). Then, when two 
vibrons are located onto the same site (m = 0), Eq.(19) 
is expressed as 

(2w - 2i fe - w)* fe (0) - V2 7fc * fe (l) = (23) 

where A k = A + 2J 2 §{2) i cos(k). 

Finally, for each wave vector k, the two-vibron dynam- 
ics reduces to a one-dimensional tight-binding problem 
on a semi-infinite lattice which exhibits two defects (see 
Fig. lb). The Schrodinger equation Eqs.(21)-(23) can 
thus be solved easily to obtain the two-vibron eigenstates 
and to determine the two-vibron energy spectrum. This 
procedure is illustrated in the following section. 

V. NUMERICAL RESULTS 

In this section, the previous formalism is applied to 
characterize the two-vibron energy spectrum of an a- 
hclix protein. However, the present theory involves a 
set of parameters which has first to be discussed. 

From the literature, the harmonic dynamics of vibrons 
in a-helices is relatively well described (see for instance 
Rcfs. [16,27]). The quantum energy for an amide-I vi- 
bration is about ujq = 1665 cm -1 and a well admitted 
value for the hopping constant is J = 7.8 cm -1 . By con- 
trast, the phonon dynamics as well as the vibron-phonon 
coupling parameter are only partially known. The mass 
M of a peptide group ranges between 1.17 10 -25 kg and 
1.91 10 -25 kg, whereas the phonon force constant W is 
expected to be about 13-19.5 Nm" 1 . As a result, the 
phonon cutoff frequency varies from 87 to 137 cm -1 . In 
the same way, a typical range for the vibron-phonon cou- 
pling term is x = 35 — 62 pN. Therefore, the vibron- 
phonon coupling parameter A ranges between 12 and 
29 cm -1 and the small polaron binding energy Eb ex- 
tends from 3 to 15 cm -1 . 

In a marked contrast with the previous parameters, 
little is known about the intramolecular anharmonicity 
of each amide-I vibration. Therefore, to determine the 
range of the anharmonic parameters, we assume that 



the amide-I vibration is equivalent to the C=0 stretch- 
ing mode of a single molecule, but with a different re- 
duced mass [i. By describing both vibrations accord- 
ing to a Morse potential, it is straightforward to show 
that wo scales as l/y/jl whereas the anharmonicity A 
scales as By comparing the harmonic frequency of 
both the amide-I vibration and the C=0 stretching mode 
(u>o = 2170 cm -1 [36]), the reduced mass of the amide-I 
mode is about 1.6-1.7 times the reduced mass of the C=0 
molecule. As a consequence, since the anharmonic pa- 
rameter for the C=0 stretching mode is A = 13.3 cm -1 
[36] , the anharmonic constant of the amide-I vibration is 
about A — 8.0 cm -1 . Note that this latter value repre- 
sents an order of magnitude of the anharmonic parameter 
since, strictly speaking, the anharmonicity of the amidc- 
I vibration depends on the details of the corresponding 
intramolecular potential. In addition, our calculations 
establish that both the cubic and quartic anharmonic 
parameters can be expressed approximately in terms of 
A by using the relation 15j^/u>o w 674 w A. 

As a result, in the following of the text, the intramolec- 
ular anharmonicity will be described by a single param- 
eter, namely the anharmonic constant A. The small po- 
laron binding energy Eb is taken as a parameter whereas 
the phonon cutoff frequency O c is fixed to 100 cm -1 . 
The temperature is set to the biological temperature, i.e. 
T = 310 K, and the hopping constant is set to J = 7.8 
cm -1 . 

The influence of the intramolecular anharmonicity on 
the two-vibron energy spectrum is shown in Figs. 2, 3 
and 4 for three typical situations. In each figure, the 
spectrum is centered onto the corrected frequency 2u> 
(Eq.(ll)) and corresponds to the two-vibron dispersion 
curves drawn in the first Brillouin zone of the lattice, i.e. 

— 7T < k < 7T. 

When E B = 8 cm" 1 and A = (Fig. 2a), the two- 
vibron energy spectrum is formed by an energy contin- 
uum associated to the TVFS. This continuum is sym- 
metrically located around 2o>o with a bandwidth equal to 
31.21 cm -1 . Below the continuum, the energy spectrum 
exhibits two bands connected to two different bound 
states. The low frequency band, located below the TVFS 
over the entire Brillouin zone, refers to bound states de- 
noted TVBS-I. The bottom of this band is located at 
24.74 cm -1 below the center of the continuum and its 
bandwidth is about 8.7 cm -1 . The second band, which 
refers to bound states denoted TVBS-II, is located be- 
low the continuum at the end of the Brillouin zone, 
only. The band disappears inside the continuum when 
I k \< k c = 2.06. When the anharmonic parameter is 
set to A = 10 cm" 1 (Fig. 2b), the TVFS bandwidth 
is reduced to 28.40 cm -1 . Note that the center of the 
continuum is rcdshiftcd due to the dependence of 2a>o 
(Eq.(ll)) on the anharmonic parameter A (not drawn in 
the figure). Then, the intramolecular anharmonicity is 
responsible for a strong redshift of the TVBS-I band as 
well as for a decrease of the corresponding bandwidth. 
The bottom of the band is located at 43.39 cm -1 below 
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the center of the continuum and the bandwidth is equal 
to 2.78 cm -1 . Finally, although the anharmonicity does 
not significantly change the shape of the TVBS-II band, 
it leads to a decrease of the critical wave vector for which 
its disappearance arises since k c = 0.70. 

When Eb is set to 10 cm -1 , the two-vibron energy 
spectrum is shown in Figs. 3a and 3b. In the harmonic 
case, i.e. A = (Fig. 3a), the energy spectrum ex- 
hibits the same features as in Fig. 2a. Nevertheless, the 
TVFS bandwidth is smaller and equal to 26.24 cm -1 . 
The width of the TVBS-I band has been reduced to 6.04 
cm -1 and the bottom of the band is located at 26.04 
cm -1 below the center of the continuum. In addition, 
the critical wave vector connected to the disappearance 
of the TVBS-II has been reduced to k c = 1.70. When 
the anharmonic parameter is set to A = 10 cm' 1 , we 
observe the same behavior as in Fig. 2b. Indeed, both 
the TVFS and the TVBS-I bandwidths decrease to reach 
23.30 cm -1 and 1.70 cm -1 , respectively. In addition, the 
bottom of the TVBS-I band is strongly redshiftcd and 
is located at 47.47 cm -1 below the center of the contin- 
uum. However, a new process arises since the intramolec- 
ular anharmonicity yields the occurrence of the TVBS-II 
band over the entire Brillouin zone. The corresponding 
bandwidth is equal to 1.53 cm -1 and the bottom of the 
band is located at 13.42 cm -1 below the center of the 
continuum. 

The situation corresponding to a strong small polaron 
binding energy is finally illustrated in Figs. 4a and 4b 
for Eb=14: cm -1 . The main difference when compared 
with the two previous cases is that the dressing effect 
is strong enough to induce the occurrence of the TVBS- 
II band over the entire Brillouin zone, even when A = 
(Fig. 4a). As previously, the intramolecular anharmonic- 
ity reduces both the TVFS continuum and the TVBS-I 
bandwidth and yields a strong redshift of the TVBS-I 
band. Moreover, it increases the frequency difference be- 
tween the TVBS-II band and the center of the continuum 
and reduces the TVBS-II bandwidth. 

The difference between the two bound states is illus- 
trated in Fig. 5. Figs. 5a and 5b show the TVBS wave 
functions with zero wave vector and which the energy 
spectrum is displayed in Figs. 3a and 3b, respectively 
(Eb = 10 cm -1 ). By contrast, Figs. 5c and 5d refer to 
the TVBS which the energy spectrum is shown in Figs. 
4a and 4b, respectively (E B = 14 cm" 1 ). The TVBS-I 
wave function is described by open circles whereas the 
TVBS-II wave function is represented by open squares. 
When A = and E B = 10 cm -1 (Fig. 5a), the TVBS-I 
wave function is maximum when to = and decreases 
as the separating distance to between the two vibrons 
increases. However, it does not exhibit a true exponen- 
tial decrease versus to. When the anharmonicity is set to 
A = 10 cm -1 (Fig. 5b), the lattice supports two bound 
states with zero wave vector (see Fig. 3b). As previously, 
the TVBS-I wave function is maximum when to = but 
the extension of the wave function has been reduced when 
compared with the harmonic case, i.e. it decreases expo- 



nentially versus to. By contrast, the TVBS-II wave func- 
tion is maximum when m = 1 and exhibits an exponen- 
tial tail as to increases. When Eb = 14 cm -1 , the lattice 
supports two bound states with zero wave vector what- 
ever the anharmonicity (Figs. 4a and 4b). When ^4 = 
(Fig. 5c), the TVBS-I wave function, which is maximum 
when to = 0, shows a significant value when to = 1 and 
does not decrease according to a true exponential. The 
TVBS-II wave function is maximum when to = 1 and 
exhibits a rather important weight when m = 0. By con- 
trast, when the anharmonicity is set to A = 10 cm -1 
(Fig. 5d), the TVBS-I and the TVBS-II wave functions 
are clearly localized onto to = and to = 1, respectively, 
both wave functions exhibiting an exponential tail versus 

TO. 

As mentioned previously, the intramolecular anhar- 
monicity is responsible for a slight decrease of the TVFS 
bandwidth. This feature is illustrated in Fig. 6 where 
the behavior of the bandwidth is shown as a function 
of A for different values of Eb- The figure exhibits two 
distinct regimes. When E B = 0, the TVFS bandwidth 
slightly increases in a linear way as the anharmonicity 
increases. The bandwidth, equal to 62.40 cm -1 when 
A = 0, is blucshifted of about 0.36 cm -1 when A = 10 
cm -1 . In a marked contrast, for finite values of Eb, the 
TVFS bandwidth is strongly reduced and decreases as A 
increases. When E B — 9 cm" 1 and A = 0, the band- 
width is equal to 28.6 cm -1 whereas it reaches the value 
25.7 cm -1 when A = 10 cm -1 , i.e. a variation of about 
2.9 cm -1 . 

In Figs. 7a and 7b, the influence of the anharmonicity 
on the TVBS-I is illustrated for different values of Eb- 
As shown in Fig. 7a, the TVBS-I binding energy Ej, 
defined as the gap between the TVBS-I with zero wave 
vector and the bottom of the TVFS continuum, decreases 
as the anharmonicity increases. When E B = 0, Ej de- 
creases from zero according to a A 2 power law for small A 
values and reaches a quasi-linear decrease for a stronger 
anharmonicity. For finite values of Eb, Ej reaches more 
rapidly the linear regime as increasing A and the slope 
of the decrease slightly increases with E B . Note that Ei 
decreases with the small polaron binding energy. In Fig. 
7b, the behavior of the TVBS-I bandwidth as function 
of A is illustrated for different values of E B . In a gen- 
eral way, the TVBS-I bandwidth decreases as the anhar- 
monicity increases. However, for a small anharmonicity, 
the TVBS-I bandwidth exhibits a rather fast decrease 
whereas it decreases more slowly for a stronger anhar- 
monicity. When Eb = 9 cm -1 , the TVBS-I bandwidth 
is equal to 7.28 cm -1 when A = and reaches 2.67 cm -1 
when A = 8 cm -1 , i.e. almost 0.3 times smaller. 

Finally, Fig. 8 displays the behavior of the critical wave 
vector k c for which the TVBS-II band disappears. In a 
general way, k c decreases as A increases and vanishes for 
a critical value of the anharmonic parameter A = A c . In 
others words, when A > A c , the TVBS-II band is located 
below the continuum. By contrast, when A < A Cl the 
TVBS-II band appears below the continuum at the end 
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of the Brillouin zone only, i.e. when | k \> k c . Note that 
A c decreases as Eb increases. Moreover, the different 
curves clearly exhibit a critical behavior since k c scales as 
k c « (A — Ac) 1 / 2 when it approaches zero. Note that, as 
shown in the inset of Fig. 8, k c shows a similar behavior 
with respect to the temperature. Indeed, for A = 8 cm -1 
and Eb — 10 cm -1 , the critical wave vector decreases as 
increasing the temperature and vanishes when T = T c = 
236 K. Therefore, when T > T c , the TVBS-II band 
is located below the continuum over the entire Brillouin 
zone. Note that k c shows the same critical behavior when 
it approaches zero and scales as k c w (T — Tc) 1 / 2 . 

VI. DISCUSSION 

To interpret and discuss the previous numerical results, 
let us first focus our attention onto the influence of the 
anharmonicity on the TVFS continuum. As shown in 
Figs. 2, 3, 4 and 6, the intramolecular anharmonicity is 
responsible for a rcdshift of the TVFS bandwidth, typ- 
ically of about a few wave numbers. Note that, when 
Eb = 0, i.e. when there is no coupling with the phonon 
bath, the anharmonicity leads to a very small blueshift 
of the continuum bandwidth. These features originate in 
the modification of the hopping constant J due to both 
the anharmonicity and the dressing effect. Indeed, as 
discussed in Sec. IV, a TVFS corresponds to an inde- 
pendent motion of the two vibrons according to an ef- 
fective hopping constant Ji$(l) = J\exp[— (1 + Ar])S(T)] 
which differs from the constant involved in the harmonic 
approximation. First, the intramolecular anharmonic- 
ity slightly increases the harmonic hopping constant (see 
Eq.(7)) yielding a blueshift of the TVFS continuum when 
Eb = 0. Then, the anharmonicity reinforces the role 
played by the coupling constant S(T) which character- 
izes the dressing effect (S(T) -> (1 + 4?y)S(T)). There- 
fore, it enhances the dressing effect and favors a decrease 
of the effective hopping constant. Although both con- 
tributions act in an opposite way, the enhancement of 
the dressing effect appears to be more efficient for non 
vanishing values of the small polaron binding energy. In 
other words, anharmonic vibrons are more sensitive to 
the dressing effect than harmonic vibrons. For instance, 
when T = 310 K, E B = 14 cm" 1 and J = 7.8 cm" 1 , the 
effective hopping constant Ji$(l) is equal to 2.32 cm -1 
within the harmonic case whereas it is redshifted to 1.96 
cm -1 when A = 10 cm -1 , i.e. a reduction of about 15%. 

The main result of the present study concerns the mod- 
ification of the bound states due to the interplay between 
the intramolecular anharmonicity and the strong vibron- 
phonon coupling. Indeed, for the harmonic lattice, the 
dressing effect is responsible for the occurrence of two 
different bound states, i.e. TVBS-I and TVBS-II. The 
TVBS-I are located below the TVFS continuum over the 
entire Brillouin zone. By contrast, for TVBS-II, two sit- 
uations occur depending on the strength of the small po- 



laron binding energy. For small values of Eb, the band 
disappears inside the continuum when the wave vector 
| k | is lower than a critical wave vector k c whereas, for 
strong values of Eb, the band is located below the contin- 
uum over the entire Brillouin zone. As when increasing 
the anharmonicity, the TVBS-I band is redshifted and 
its bandwidth is strongly reduced. In the same way, the 
anharmonicity modifies the nature of the TVBS-II band. 
If the band is resonant with the continuum, the anhar- 
monicity yields a decrease of the critical wave vector. 
Therefore, if it is strong enough, the anharmonicity al- 
lows for the TVBS-II band to get out of the continuum 
over the entire Brillouin zone. When the TVBS-II band 
is not resonant with the continuum, the anharmonicity 
induces a redshifted of the band as well as a decrease 
of its bandwidth. All the previous features are accompa- 
nied by a modification of the wave functions of the bound 
states. 

To understand these features, we can take advantage 
of the equivalence between the two-vibron dynamics and 
the tight-binding problem on the one-dimensional lattice 
displayed in Fig. lb. Within this equivalence, bound 
states in the real molecular lattice are described in terms 
of localized states in the equivalent lattice, the localiza- 
tion occurring due to the presence of the two defects. 
For a given k value, the previous results clearly show 
that the system supports one or two bound states, de- 
pending on the values taken by the relevant parameters 
of the problem (anharmonicity, small polaron binding en- 
ergy, temperature ...). In other words, if these parame- 
ters are allowed to vary, a localization transition arises in 
the equivalent lattice. Such a transition discriminates be- 
tween two "phases". The first phase corresponds to the 
presence of a single localized state (i.e. a single bound 
state in the real lattice) whereas the second phase is con- 
nected to the occurrence of two localized states (i.e. two 
bound states). A localized state is characterized by its lo- 
calization length £ which refers to the length of the bond 
between two vibrons in a bound state. The disappear- 
ance of a localized state is accompanied by a divergence 
of the localization length since it now refers to two inde- 
pendent vibrons. 

Such a process is similar to a critical transition and 
consequently, we can use the tools of the Renormaliza- 
tion Group (RG) theory to understand the transition and 
predict the occurrence of localized states (i.e. of bound 
states) . This can be achieved by performing a decimation 
of the Schrodinger equation as detailed in Refs. [37,38] 
and summarized in appendix B. The decimation consists 
in eliminating from the initial equivalent lattice one site 
over two to arrive at a scaled lattice with twice the lat- 
tice spacing. If the initial lattice is at a critical point, i.e. 
if the localization transition takes place, the localization 
length £ is infinite and the dynamics is in a self-similar 
situation. As a result, no change in the critical parame- 
ters accompanies the length scaling and the scaled lattice 
remains at a critical point. The decimation procedure 
can be applied recursively until the lattice parameter has 
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been increased up to infinity. By operating in such a way, 
we decrease drastically the number of sites in the equiv- 
alent lattice, and obtain a scaled lattice which reduces to 
the two side sites m = and m — 1, only (see Fig. lb). 
It is thus possible to study exactly the Schrodinger equa- 
tion of such a critical system and then to characterize 
the critical properties of the initial lattice. By following 
this procedure (see appendix B), the condition for the oc- 
currence of localized states in the equivalent lattice, i.e. 
for the occurrence of bound states in the real lattice, is 
defined as 

(A k -T k )(B-T k ) = 7fe 2 (24) 

where the parameters occurring in Eq.(24) are defined in 
Sec. IV. 

As shown in Fig. 9, Eq.(24) defines the critical curve 
in the space of the parameters which separate the phase 
I, with a single bound state, from the phase II, with two 
bound states. In phase I, the nature of the bound state 
depends on the relative values of the parameters and 
three main situations occur. When Ak » B, the lo- 
calization arises around the first site m = 0. In other 
words, TVBS-I characterizes the trapping of two vibrons 
around the same amidc-I vibration. By contrast, when 
Ah « B, the localization occurs on the second site 
m = 1 and the two vibrons are trapped onto nearest 
neighbor amide-I modes. In the intermediate case, when 
Ak ~ B, the localized state is a superimposition of states 
localized on sites m = and m — 1. As a result, TVBS- 
I characterizes a superimposition of states trapped over 
both the same amide-I mode and two nearest neighbor 
amidc-I modes. In phase II, the lattice supports two 
bound states. When Ak » B, TVBS-I corresponds to 
the trapping of the two vibrons onto the same amidc- 
I mode whereas TVBS-II characterizes the trapping of 
the two vibrons onto nearest neighbor amide-I modes. 
As shown in Fig. 9, the revert situation takes place 
when A k « B. Finally, when A k < B < 2A k , both 
bound states appear as combinations of states involving 
the trapping of the two vibrons onto the same amide-I 
mode and onto nearest neighbor amide-I modes. 

In the harmonic a-helix protein, i.e. A — 0, Eq.(ll) 
yields Ak = B = Eb- Therefore, the critical curve 
Eq.(24) shows that when E B < 4J$(0)cos(fc/2), the lat- 
tice exhibits a single bound state which mixes the trap- 
ping onto the same amide-I mode and onto two nearest 
neighbor amide-I modes (see Fig. 5a). By contrast, when 
Eb > 4J$(0)cos(fc/2), the lattice supports two bound 
states, both being a superimposition of states involving 
the trapping onto the same and onto two nearest neigh- 
bor amide-I modes (see Fig. 5c) . Note that the condition 
for the occurrence of two bound states over the entire 
Brillouin zone is Eb > 4J$(0). As when increasing the 
anharmonicity A, Eq.(ll) shows that Ak > B. For in- 
stance, when T = 310 K, E B = 14 cnT 1 and A = 10 
cm -1 , A is equal to 28.29 cm" 1 whereas B is equal to 
16.54 cm -1 . Therefore, the anharmonicity decreases the 



hybridization between the two kinds of trapping. As a 
result, TVBS-I corresponds essentially to the trapping 
onto the same amide-I mode whereas TVBS-II, when it 
is present, characterizes the trapping onto two nearest 
neighbor amidc-I modes (see Figs. 5b and 5d). Note that 
the hybridization between the two trapping processes is 
also reduced when the hopping constants Tk and jk de- 
crease, i.e. when the different sites m become uncorre- 
cted. This feature arises when the TVBS wave vector 
approaches ir as well as when the temperature increases 
due to the dressing effect. 

Although the phase diagram displayed in Fig. 9 al- 
lows for a complete understanding of the nature of the 
bound states, it involves parameters which are not in- 
dependent. For instance, both Ak and B depend on A 
and Eb- However, from Eq.(24), we can define a phase 
diagram in the (A, Eb) parameter space. Such diagram 
is illustrated in Fig. 10 for a zero wave vector. The tem- 
perature is fixed to T = 310 K whereas three different 
values for the hopping constant J have been considered. 
The curve discriminates between a phase with a single 
bound state and a phase with two bound states, both 
states being located below the continuum over the entire 
Brillouin zone since k — 0. The figure clearly shows that 
the anharmonicity favors the occurrence of two bound 
states by decreasing the value of the required small po- 
laron binding energy. For instance, when J = 7.8 cm -1 
and A — 0, E B must exceed 11.5 cm -1 to allow for the 
occurrence of the TVBS-II. By contrast, this value is re- 
duced to 8.7 cm -1 when A = 8 cm -1 . Note that the 
curve is pushed down as decreasing the hopping constant 
as well as increasing the temperature due to the dressing 
effect. 

From a physical point of view, the previous phase di- 
agrams can be easily understood in terms of localization 
in the equivalent lattice. Indeed, since the equivalent lat- 
tice exhibits two defects onto the sites m — and m = 1, 
it can supports two states localized onto the sites m = 
and m = 1. However, due to the hopping processes, these 
two localized states overlapp and interact to generate new 
localized states which appear as superimpositions of the 
two previous states. Such a mechanism is accompanied 
by a splitting of the frequency of the new localized states 
which depends on both the strength of the coupling be- 
tween the original states as well as on their frequency. 
As a result, we obtain a high frequency localized state 
(TVBS-II) and a low frequency localized state (TVBS- 
I). If the splitting is strong enough, the high frequency 
localized state is pushed back into the continuum and a 
single localized state remains. By contrast, if the split- 
ting is rather weak, both localized states can be located 
below the continuum. As shown in Eqs.(21)-(22), the in- 
tramolecular anharmonicity increases the energy differ- 
ence between the two original localized states. In addi- 
tion, since it enhances the dressing effect, it reduces the 
coupling between these two states. As a result, the an- 
harmonicity favors a weak hybridization between the two 
original localized states. The main consequence is that 
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TVBS-I and TVBS-II are essentially localized around 
m = and m = 1, respectively. 

The present paper has clearly established the inter- 
play between the intramolecular anharmonicity and the 
strong vibron-phonon coupling. We have shown that the 
anharmonicity modifies the vibron-phonon interaction 
which results in an enhancement of the dressing effect. 
Therefore, anharmonic vibrons are more sensitive to the 
dressing than harmonic vibrons and are characterized by 
smaller effective hopping constants. In addition, we have 
shown that both nonlinear sources break the vibron inde- 
pendence and favor the occurrence of two kinds of bound 
states which the properties strongly depend on the inter- 
play between the anharmonicity and the dressing effect. 
Such a dependence was summarized in a phase diagram 
which characterizes the number as well as the nature of 
the bound states as a function of the values taken by the 
relevant parameters of the problem. In the harmonic sit- 
uation, the two bound states appear as combinations of 
states involving the trapping of the two vibrons onto the 
same amide-I mode and onto nearest neighbor amidc- 
I modes. By contrast, the intramolecular anharmonic- 
ity reduces the hybridization between these two kinds of 
trapping so that the low frequency bound state refers to 
the trapping of the two vibrons onto the same amide-I 
mode whereas the high frequency bound state character- 
izes the trapping onto nearest neighbor amide-I vibra- 
tions. In addition, the anharmonicity strongly reduces 
the dispersion of the bound states and thus enhances 
their breather-like behavior [19]. 

To conclude, let us mention that forthcoming works 
will be devoted to the fundamental question of the TVBS 
lifetime due to the coupling AH (Eq. (17)) with the re- 
maining phonons. Such a problem was studied in a recent 
paper [18] in which the decay of TVBS into either bound 
or free states was described by considering a weak vibron- 
phonon coupling. However, in the present context, both 
the anharmonicity and the dressing effect modify the na- 
ture and the number of the bound states as well as the 
coupling with the remaining phonons. As a consequence, 
a new theory must be performed to characterize the dif- 
ferent pathways for the decay of the TVBS. 

Appendix A: Unitary transformation to remove 
the intramolecular anharmonicity 

By assuming that the cubic anharmonicity is one order 
of magnitude greater than the quartic anharmonicity, a 
perturbative parameter A is introduced so that the vibra- 
tional Hamiltonian of the nth amide-I mode is written as 



h n = h.M +\V n + \ 2 W n + .... 



(Al) 



where = u Q b+b„, V n = 73(6^ + 6 n ) 3 and W n = 
74(6+ + b n ) 4 . The anharmonic terms in Eq.(Al) can be 
removed by performing a perturbative unitary transfor- 
mation T n = exp(S n ), where S n is expanded as 



Under this transformation, the Hamiltonian h n becomes 
h n = h^+X(V + [Si 1 \h^}) + 
+X 2 (W+[S^\h^} + [s£\V\ + [SW.fcW]]) 
+... " (A3) 

Since the required unitary transformation must diagonal- 
ize the Hamiltonian, the diagonal terms in Eq.(A3) par- 
ticipate in the diagonalization scheme of the Hamiltonian 
whereas the non diagonal terms, which must vanish, al- 
low us to determine the S n operator order by order. At 
the second order with respect to the anharmonic param- 
eters, we obtain 

SP = 11 K 3 -b 3 n + 9(bt 2 b n b+bl + b+ b n )] (A4) 

W 

and 

s^-\[^ + ^?][b^-bi] + 

{ 1± _ 3 (^) 2 ][26+ 3 &„ _ 2b + b 3 + 35 +2 _ 36 2j (A5) 
Wo W 

which lead to the renormalizcd Hamiltonian 



(w - 2A)blK 



Ab+b+b n b n 



(A6) 



where the irrelevant constant was disregarded and where 
A denotes the positive anharmonic parameter defined in 
Eq.(6). 

At this step, the full vibron-phonon Hamiltonian H can 
be expressed in an improved way by applying the general 
unitary transformation T = J\ n T n . The transformation 
modifies the lateral interaction as well as the vibron- 
phonon coupling Hamiltonian leading to the occurrence 
of vibron-conserving terms and vibron-nonconserving 
terms. Nevertheless, since the internal frequency w is 
more than two orders of magnitude greater than the an- 
harmonic parameters 73 and 74, the nonconserving terms 
are weak when compared with the conserving terms and 
will be neglected. Such a procedure requires the knowl- 
edge of the transformation of the vibron displacement 
operator b n + 6+ expressed as 

b+ + b n = b+ + b n + 2^ [b+ 2 +bl- 6b+b n - 3] 



w 



+ [22(^)2 _ Qj± ][b + + bn + b +2 bn + b + b 2 n] 

ujo wo 

+ [3(^)2 + ^+3^3] 
WO Wo 



.74- 



(A7) 



By using Eq.(A7), the transformation of the operators 
(b+ + b n ) 2 and (6+ + b n )(b+ +1 + 6 n +i) involved in the 
vibron-phonon Hamiltonian can be determined easily. 
By restricting the calculations to the vibron-conserving 
terms, we thus obtain 



(b+ + b n f = 1 + 88(^) 2 - 12^ + 2(1 + 2 V )btb n 
wo w 



,74 



5„ = A5( 1 )+A 2 5( i 2 ) + ... 



(A2) 



+2r]b^b+b n b n + nonconserving terms 



(A8) 
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where 77 is defined in Eq.(7), and 

(bt + b n )(b+ + 6„0 = 72(^) 2 [6+6„ + b+b n ,] 
+ [1 + Ui llf-12l±]{b+b n ,+b+b n } 

Wo ^0 

+ [22(^r-122±]b+(b+b n + b+b n ')b nl 
+ [22(^-12^]b+(bib n + b+b n ,)b n 

UJ UJQ 



+^?K 2 bl, + b+?bl] + iu{J±Y[b+b+,b n b n ,] 



+nonconscrving terms 



(A9) 



Finally, by using Eqs.(A7)-(A9), the tranformed 
Hamiltonian Eq.(5) is obtained straightforwardly. 

Appendix B: Decimation of the two-vibron 
Schrodinger equation 

For each k value, the two-vibron Schrodinger equation 
Eqs. (21)-(23) can be reduced by using a decimation pro- 
cedure [37,38]. To proceed, let us rewrite the Schrodinger 
equation Eqs. (21)-(23) as 

(A + 2a)* fe (0) = -V2g* k {l) (Bl) 

(A + 6)* fc (l) = -V2g^ k (0) - * fc (2) 

A* fe (m) = -* fc (m+ 1) - * fe (m- l);m= 2,3,4,... 

where A = (w — 2u) a )/T k , a = A k /T k , b = B/T k and 
9 = 7k AV 

In Eq.(Bl), we eliminate the wave functions connected 
to the even sites by substituting their expressions in the 
Schrodinger equation of the odd sites. Nevertheless we 
keep unchanged the first two Schrodinger equation asso- 
ciated to the defect sites m — and m = 1. We thus 
obtain a new set of Schrodinger equations for the odd 
sites only, as 

(2 — A 2 — 2(Aa + l))* fc (0) = V2\g* k (l) (B2) 

(2 - A 2 - (Xb + l))* fc (l) = -V2A. 9 * fc (0) - * fc (3) 

(2 - A 2 )* fe (m) = -* fc (m + 1) - * fc (m - 1); m = 3, 5, 7, ... 

Eq. (B2) characterizes the Schrodinger equation of the 
rescaled lattice with new parameters defining the RG 
transformation 



A« 



2 -A 2 
a y ~' = — Xa — 1 
6« = -A6-l 

9 {1) = -x 9 



(B3) 



The previous decimation procedure allows us to de- 
fine the critical values of the parameters responsible for 
the occurrence of bound states. Indeed, bound states 
correspond to a localization of the wave function ty k (m) 
around the sites m — and m = 1. Therefore, when the 



decimation is applied recursively p times, the neighboring 
site of the site m = 1 is pushed to infinity and we thus 
obtain an ultimate scaled lattice formed by the two side 
sites m = and m = 1, only. The Schrodinger equation 
can thus be solved exactly and its two eigenvalues are 
expressed as 



A (p) = _ 



2a( p ) + M p ) 



±\ ( 



2 a (p) + &(p) 



) 2 - 2(o(p)&(p) - g^ 2 ) (B4) 



A TVBS occurs when its frequency is at least equal to 
the minimum of the TVFS continuum. This condition is 
obtained when A = A c = — 2. At the critical point A = A c , 
the scaled values of the parameter A satisfy A^ = ... = 
\(p) = \ c _ This parameter A becomes scale invariant 
and is a fixed point of the RG transformation. Indeed, 
when the initial lattice dynamics is at a critical point, the 
localization length £ is infinite. No change in the critical 
parameters accompanies the length scaling so that the 
scaled lattice remains at a critical point. 

From Eq.(B3), the scaled values of the parameters at 
the critical point are expressed as 



X (p) 

n(P) 



2 p (a- 1) + 1 
=2 p (fo-l) + l 

9 { 



» = 2 p g 



(B5) 



Combining Eqs.(B5) and (B4) for p tends to infinity leads 
to the conditions for the occurrence of a bound states in 
terms of the reduced parameters a,b and g as 



(a-l)(b-l)=g 2 



(B6) 



Eq.(B6) is equivalent to Eq.(24) with the corresponding 
definition of a, b and g. 
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Figure Caption 

Figure 1 : (a) Equivalence between the two-vibron 
Schrodinger equation and the dynamics of a single fic- 
titious particle moving quantum mechanically on a 2D 
lattice, (b) For each wave vector k, the two-vibron dy- 
namics reduces to a one-dimensional tight-binding prob- 
lem on a semi-infinite lattice which exhibits two defects 
(see the text). 



Figure 2 : Two-vibron energy spectrum for Eb = 8 
cm" 1 , J = 8 cm" 1 , Q c = 100 cm" 1 , T = 310 K and for 
(a) A = cmT 1 and (b) A = 10 cm" 1 . 

Figure 3 : Two-vibron energy spectrum for Eb = 10 
cm" 1 , J = 8 cm- 1 , Q c = 100 cm" 1 , T = 310 K and for 
(a) A = enr 1 and (b) A = 10 cm" 1 . 

Figure 4 : Two-vibron energy spectrum for Eb = 14 
cm^ 1 , J = 8 cm" 1 , Q c = 100 cm" 1 , T = 310 K and for 
(a) A = cm -1 and (b) A = 10 cm- 1 . 

Figure 5 : Zero wave vector two-vibron bound state 
wave functions for J = 8 cm -1 , O c = 100 cm -1 , T = 310 
K, and (a) E B = 10 cm" 1 , A = cm" 1 , (b) E B = 10 
cm" 1 , A = 10 cm -1 , (c) E B = 14 cm -1 , A = cm -1 , (d) 
Eb = 14 cm -1 , A = 10 cm -1 . Open circles represent the 
TVBS-I wave function whereas open squares characterize 
the TVBS-II wave function (see the text). 

Figure 6 : Behavior of the two-vibron free state band- 
width as a function of the anharmonicity A for J = 8 
cm- 1 , n c = 100 cm" 1 , T = 310 K and E B = cm" 1 
(full circles), E B — 5 cm -1 (full squares), E B = 7 cm -1 
(full triangles), Eb — 9 cm -1 (open circles), and Eb = 11 
cm -1 (open squares). 

Figure 7 : (a) TVBS-I binding energy and (b) band- 
width as a function of the anharmonicity A for Eb = 
cm- 1 (full circles), E B = 5 cm -1 (full squares), E B = 7 
cm -1 (full triangles), E B = 9 cm -1 (open circles), 
Eb = 11 cm- 1 (open squares) and for J = 8 cm -1 , 
n c = 100 cm" 1 and T = 310 K. 

Figure 8 : Behavior of the critical wave vector for which 
the disappearance of the TVBS-II band takes place as 
a function of the anharmonicity A for E B = 8 cm -1 
(full line), E B = 9 cm" 1 (dotted line), E B = 10 cm" 1 
(short dashed line), Eb = 11 cm -1 (medium dashed line) 
(J = 8 cm" 1 , n c = 100 cm" 1 , T = 310 K). The inset 
represents the evolution of the critical wave vector as a 
function of the temperature. 

Figure 9 : Phase diagram in the parameter space. The 
critical curve discriminates a phase with a single bound 
state from a phase with two bound states (see the text). 

Figure 10 : Phase diagram in the (A, E B ) parameter 
space for a zero wave vector (fi c = 100 cm -1 , T = 310 
K). The critical curve discriminates a phase with a single 
bound state from a phase with two bound states (see the 
text). 
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